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Abstract 

We describe our recent progress in modeling supernova-fallback disks and neutron 
star (NS) atmospheres. 

We present a first detailed spectrum synthesis calculation of a SN-fallback disk 
composed of iron. We assume a thin disk with a radial structure described by the 
a-disk model. The vertical structure and emission spectrum are computed self- 
consistently by solving the structure equations simultaneously with the radiation 
transfer equations under non-LTE conditions. We describe the properties of a spe- 
cific disk model and discuss various effects on the emergent UV/optical spectrum. 

We investigate Compton scattering effects on the thermal spectrum of NSs. In 
addition, we constructed a new generation of metal line-blanketed non-LTE model 
atmospheres for NSs. It is compared to X-ray burst spectra of EXO 0748— 676. It 
is possible that the gravitational redshift, deduced from absorption lines, is lower 
(z=0.24) than hitherto assumed (z=0.35). Accordingly, this would result in NS mass 
and radius lower limits of M >1.63 M Q and i?>13.8 km. 
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1 Introduction 



We report on our recent progress in radiation transfer modeling for supernova- 
fallback disks and non-magnetic neutron star (NS) atmospheres. Disk model- 
ing is described in Sect. 2. We then present effects of Compton scattering on 
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the thermal spectra of NS atmospheres (Sect. 3). In Sect. 4 we describe the 
construction of new line-blanketed NLTE atmosphere models for NSs. These 
models are applied to X-ray burst spectra of EXO 0748— 676 (Sect. 5). 



2 NLTE models for SN fallback disks 



Anomalous X-ray pulsars (AXPs) are slowly rotating young isolated NSs. 
Their X-ray luminosities greatly exceed the rates of rotational energy loss. 
It is now generally believed that AXPs are magnetars and that their X-ray 
luminosity is powered by magnetic energy. As an alternative explanation the 
X-ray emission was attribut ed to accretion from a disk that is made up of 
supernova-fallback material (le.g. van Paradijs et all Il995l ). It is still a matter 
of debate whether the disk model is appropriate. For example, the discovery 
of optic al pulsations in 4U 014 2+61 with the same period like the X-ray pul- 
sations (IKern fc Martini . |2002| ) appears to be a strong argument against the 
disk model. It was argued that reprocessing of the pulsed NS X-ray emission 
in the disk cannot explain the high optical pulsed fraction, because disk radi- 
ation woul d be dominated by vis cous dissipati on and not by reproce ssed NS 
irradiation (IKern fc Martini . 120021 ) . In contrast. lErtan fc Chengl (120041 ) showed 
that the optical pulsations can be explained either by the magnetar outer gap 
model or by the disk-star dynamo mod el. A spectral break in t he optical spec- 
trum of 4U 0142+61 was discovered by iHulleman et al.l (120041 ) and also taken 
as an argument agai nst the disk model. The recent discovery of mid-IR emis- 
sion from this AXP (jWang et al.l . 120061 ) . however, has strongly rekindled the 



interest in studies of fallback-disk emission properties. The emis sion was at 



tribut ed to a cool, passive dust debris disk, but it was shown by lErtan et al 
( 120061 ) that it can be explained by an active, dissipating gas disk. 



The emission from fallback disks was hitherto modeled with blackbody spec- 
tra. In view of the importance of disk models for the quantitative interpre- 
tation of observational data it is highly desirable to construct more realistic 
models by d etailed radiation-t ransfer calculations. We employ our computer 
code AcDc ( jNagel et all 120041 ). For a detailed descrip tion of the calculation 
of fallback disk models we refer to lWerner et al.l ( 20061 ) . The radial disk struc- 
ture is calculated assuming a stationary, Keplerian, geometrically thin a-disk. 
For the results reported here we selected the following model parameters. The 
NS mass is 1.4 M . The radii of the inner and outer disk edges are 2000 and 
200 000 km, respectively. The accretion rate was set to M = 3 • 10~ 9 M /yr, 
which is a relatively large value. It is motivated by the fact, that this cor- 
responds to the upper limit for the non-detection of UV radiation from a 
possible disk a round SN 1987A un der the assumption of an outer disk radius 
of 100 000 km ([Graves et all 120051 ). Our code allows for the irradiation of the 
disk by the central source, however, the results presented here are computed 
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Fig. 1. Spectrum of a fallback-disk model seen under three inclination angles i. 
Broad iron-line blends are detectable even in the almost edge-on case. 



with zero incident intensity. The disk is represented by nine rings. The verti- 
cal structure of each ring is determined from the simultaneous solution of the 
radiation transfer equations plus the structure equations (radiative and hydro- 
static equilibrium). They also consist of the non-LTE rate equations for the 
atomic population densities. The chemical composition of SN fallback material 
is not exactly known. It depends on the amount of mass that goes into the 
disk. A disk with a small mass (s ay < 0.001 M ) will be composed of silicon- 
burning ash (IMenou et all l200ll ). For simplicity we assume a pure-iron disk 
composition but also tested a composition that represents silicon-burning ash. 
It contains iron (80% mass fraction) as well as silicon and sulfur (10% each). 
Fig. 1 displays the overall disk spectrum. We summarize the model properties 
as follows: 



• Depending on the inclination, the disk flux can be a factor of two higher or 
lower compared to a blackbody radiating disk at the same inclination. 

• Strong iron line blanketing causes broad (>100 A) spectral features that 
could be detectable even from almost edge-on disks. 

• Limb darkening affects the overall disk spectrum (in addition to the geo- 
metric foreshortening factor). Depending on inclination and spectral band, 
the disk intensity varies up to a factor of three. 

• The overall disk spectrum is independent of the detailed chemical com- 
position as long as Fe is the dominant species. In particular, a pure-Fe 
composition is spectroscopically indistinguishable from Si-burning ash. 

• The overall disk spectrum is hardly influenced by non-LTE effects, however, 
equivalent widths of individual line blends can change by a factor of two. 
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Fig. 2. Emergent spectra of three pure-H models with different T e g- (1, 3, and 5-10 6 K) 
computed with Thomson (dashed lines) and Compton (full lines) electron scattering. 



3 Compton scattering in NS model atmospheres 



The nature of isolated NS (INS) surface layers is not exactly known. At some 
conditions (depending on surface temperature, magnetic field strength and 
chemi cal composition) a surface can b e solid, liquid, or have a plasma enve- 
lope (lLai fc Salpeterl . Il997l : Q l200ll ). In the last case the envelope can be 



considered as a NS atmosphere, and the structure and emergent spectrum of 
this atmosphere can be computed by using stellar model atmosphere meth- 
ods. Such modeling was performed by many groups and the model spectra 
were used to fit observed INS X-ray spectra. One of the important results of 
these works is that spectra of light elements (hydrogen and helium) atmo- 
spheres with low magnetic field are significantly harder than the correspond- 
ing blackbody spectra. These elements are fully ionised in atmospheres with 
T c s > 10 6 K. Therefore, the true opacity (mainly due to free-free transitions) 
decreases with photon energy as E~ 3 . At high energies electron scattering 
is larger than the true opacity and photons emitted deep in the atmosphere 
(where T > T e fj) escape after few scatterings on electrons. In all previous work 
concerning INS, coherent (Thomson) electron scattering is considered. As a 
result, emergent spectra are very hard. But such a situation is very favorable 
to change the photon energy due to Compton down-scattering. 

It is well known that Compton down-scattering determines the shape of emer- 
gent model spectra o f hotter NS atmo s phere s with T cS ~ 2 • 10 7 K close to 



the Eddington limit (le.g. Pavlov et all Il99ll ). These model spectra describe 



the observed X-ray spectra of X-ray bursting NS in low-mass X-ray binaries 
(LMXBs), and they are close to diluted blackbody spectra with a hardness 
factor f c ~ 1.5-1.9. But models with Compton scattering taken into account 
were not calculated for relatively cool atmospheres with T e s < 10 7 K. There- 
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fore, at present time, the effect of Compton scattering on emergent spectra of 
INS model atmospheres with T e g- < 5 • 10 6 K is not well known. 

We investigated the effect of Compton scattering on hot (T eS > 10 6 K) INS 
atmospheres with weak magnetic fields (B < 10 8 G). In order to compute 
LTE model atmospheres in hydrostatic and radiative equilibrium we solve 
the radiation transfer equation with the Kompaneets operator. We calcu- 
lated a set of models with T e g in the range 1-5- 10 6 K, with two values of 
surface gravity (logq=13.9 a nd 14.3), and different chemical compositions 



(jSuleimanov &: Werner! . 120071 ) and the results can be summarized as follows: 



Radiation spectra computed with Compton scattering are softer than those 
computed without Compton scattering at high energies (E > 5 keV) for 
light element (H or He) model atmospheres (Fig. 2). 

The Compton effect is most significant in H model atmospheres and models 
with low surface gravity. 

Compton scattering is less important in models with solar abundance of 
heavy elements. 

The emergent spectra of the hottest (T e g > 3 • 10 6 K) models can be de- 
scribed by diluted blackbody spectra with hardness factors ~ 1.6-1.9. 



4 Metal line blanketed NLTE models for NS atmospheres 



We are computing a new grid of model atmospheres for weakly magnetized 
INSs. Since it covers a large T ef j-range up to high temperat ures (1-10-10 6 K) it 



takes into account deviations from LTE for the first time. I Werner fc Deetjen 



(120001 ) have shown that NLTE effects can affect the temperature structure of 
the outer photospheric layers. The models are line-blanketed and calculated for 
different chemical compositions comprising H, He, C, N, O, and iron (Rauch, 
Suleimanov, Werner, in pre p.). We employ the TMAP (Tubingen Model At- 



mosphere Package) codes (IWerner et all 12003!) and account for opacities of 



millions of iron lines in a statistical way (IRauch fc Deetjenl . 120031 ) . The model 
atoms are constructed using energy levels from NIST and oscillator strengths 
and photoionisation cross-sections calculated by the Opacity Project (TIP- 
TOPbase). As an example, in Fig. 3 we display the X-ray flux spectra of 
models with different T e g, solar element composition, and log g=14.39. 



5 Model comparison to X-ray burst spectra of EXO 0748— 676 



Cottam et al.l (120021 ) have announced the detection of absorption features in 



the burst spectra of the LMXB EXO 0748—676. The three most significant 
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Fig. 3. Spectra of solar composition NLTE model atmospheres with increasing T c g 
(1-10-10 6 K in steps of 10 6 K) at log 5=14.39. 

features were identified as photospheric absorption lines, redshifted by z=0. 35: 
A feature at 13.00 A in the early-burst phases was assigned to the Fexxvi 
n=2-3 transition at An=9.5-9.7 A. In the late-burst phase a feature appears 
at 13.75 A and it was assigned to the Fexxv n=2-3 transition at Ao=10.2 A. 
A double feature at 25.2/26.0 A, also appearing in the late-burst phase, was 
interpreted as a broad self-reversed line profile from O vin Ly a (Aq=18.97 A). 



We can use our models to check these line identifications. The upper curve in 
Fig. 4 is the spectrum of a Fe-dominated (99%) model with T eff =8 ■ 10 7 K, 
logq=14 .39, redshifted by z= 0.35. This Figure can be directly compared to 
Fig. 1 in lCottam et al. fj2002h . It is obvious that the aforementioned Fexxvi/ 
Fe xxvii lines at (redshifted) A=13.0/13.75 A are much too weak in the model. 
In models with different T e ff t hese lines are even weaker, so that the identifi- 
cation by Cottam et al.l ( 20021 ) is not confirmed. 



Other Fe lines in the models are much stronger and we suggest another possible 
line identification in the observed spectra. The lower curve in Fig. 4 shows a 
solar-composition model with the same T e g and logg but with a different 
redshift, namely z=0.24. The observed lines at A= 13.0/13. 75 A could stem 
from Fexxiv (n=2-3, A =10.6-11.4 A). The O vm Ly a line is located at 23.6 A 
in the 2=0.24 redshifted model and a weak feature is seen in the observed 
spectrum at this location. We would like to stress that our analysis is in a 
very preliminary stage. But our models clearly indicate that a re-analysis of 
the X-ray burst spectra of EXO 0748—676 is necessary. 
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Fig. 4. Spectra from two models with z=0.35 (top) and z=0.24 (bottom). They 
suggest that the observed absorption lines in EXO 0748— 676 at A=13.0/13.75 A 
could stem from Fexxiv at z=0.24 rather than from Fexxvi/ Fexxvn at z=0.35. 
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Fig. 5. Allowed values for M and R of EXO 0748-676 for redshifts z=0.24 
and z=0.35 (straight lines; thick portions of the graphs d enote the mass rang e 
1.4-1.8 Mq) compared to various theoretical M-R relations ( Haensel et al. . 2006). 
The t hick dot on the z=0.35 line denotes the minimum M and R derived by Ozell 
(2006). The arrow indicates the shift of this result when we assume z=0.2A. 



Fig. 5 shows the allowed values for NS mass and radius for redshifts ,2=0.24 and 
z=0.35 compared to various theoretical mass-radius relations. While z=0.35 
gives radii of R=9— 12 km for a mass-range of M= 1.4— 1.8 M Q , our red- 
shift 2=0.24 gives larger radii, namely 12-15 km, which corresponds to stiff 
equations-of-state and excludes mass-radius relations based on exotic matter. 
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This result is in line with a recent study of EXO 0748— 676 by lOzell (120061 ) 
using additional observational constraints. With z=0.35 she derives minimum 
values of mass and radius, M >2.10 M and R >13.8 km. A reduction of the 
redshift to 2=0.24 has a negligible effect on her radius determination but the 
lower mass limit is reduced to 1.63 M Q . 
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